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Abstract 

The inhibition of the total HI n «-> 1 transition rate by delayed resonant reabsorption of 
HI (n + 1) — > 1 photons by HI n — > 1 line which is possible due to cosmological redshift 
is considered semi-analytically. The method taking into account this effect in the frame 
of simple three-level approximation model of recombination is suggested. It is confirmed 
that the resonant feedbacks affect ionization fraction at the level about 0.2% for the epoch 
of last scattering. 

Similar consideration of Hel 2 l P <-> l 1 ^ =>■ Hel 2 3 P <-> l 1 ^ feedback for helium is 
provided. It is confirmed that allowance of this feedback leads to increase of predicted free 
electron fraction by 0.12% at z ~ 2300. It is shown that taking into account absorption 
and thermalization of Hel 2 X P <-> l^S resonant superequilibrium photons (during their 
redshifting to the Hel 2 3 P <-> l 1 ^ frequency) by small amount of neutral hydrogen (10~ 7 — 
10 -4 of total number of hydrogen atoms and ions) existing in helium recombination epoch 
is important for correct consideration of this helium feedback. 

Keywords: cosmological recombination, CMB, anisotropy, feedback, hydrogen, helium, 
deuterium 



1 Introduction 



Cosmological recombination plays an important role for the observations of the 
CMBR anisotropy, because the visibility of epoch at fixed z directly depends on concen- 
tration of free electrons. Today it is clear that for correct analysis of experimental data 
from Planck mission (launched at May 2009) and other future experiments on anisotropy 
measurements the ionization fraction as a function of redshift should be calculated with 
relative accuracy within 0.1% at least. Also the codes for such calculations should be fast 
enough for effective determination of cosmological parameters from CMBR anisotropy 
data. Several approaches are suggested to satisfy these conditions: three-level approxi- 
mation (TLA) with fudge factors (e.g. recfast, Seager et al. 1999, Wong et al. 2008), 
training-conception (e.g. RICO, Fendt et al. 2009), multy-level codes with not very 
large principal quantum number (5 < n max < 10), and others. In a number of works 
the different approach was described (e.g. Dubrovich and Grachev 2005) and used (e.g. 
Kholupenko and Ivanchik 2006): the TLA-model has been successfully modified to in- 
clude transitions from high excited states, helium recombination through ortho-channel, 
HI Lya =^ HI 2s<->ls feedback, and hydrogen continuum absorption. These works show 
that the way of modification of TLA-model appears fruitful and in principal in this way 
TLA-model can be developed to the form taking into account cumulative action of "thin" 
effects considered until now for precise calculations of cosmological recombination. 

The radiative feedbacks for the cosmological recombination problem have been 
considered in a number of works (Chluba and Sunyaev 2007, 2009; Switzer and Hirata 
2008). The method presented in this work is similar to the method developed by Chluba 
and Sunyaev (2009), but in present work some additional analytical approximations have 
been obtained. Use of these approximations allows us to speed up the recombination 
calculations. The main purpose of this paper is to show how the effect of HI Ly(n + 1) =>- 
HI n<^l feedbacks for hydrogen (Chluba and Sunyaev 2007, 2009) and Hel 2 1 P <-> l 1 ^ =>- 
Hel 2 3 P <-> l 1 ^ feedback for helium (Switzer and Hirata 2008, Chluba and Sunyaev 2009) 
can be taken into account in the frame of TLA-model and to provide an independent 
calculation of these effects for comparison. 

2 Cosmological model 

All calculations have been performed in the frame of standard cosmological model. Cor- 
responding values of cosmological parameters are indicated in Tab. 1. Dependence of 
Hubble constant on redshift is given by the following 

H(Z) = tfoV^A + On(l + Zf + a 6 l(l + ZY (1) 

Total concentration of atoms and ions for the certain component of plasma depends on 
redshift by power law: N ~ (1 + z) 3 . 
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Table 1: Parameters of the standard cosmological model 



Value description 



Symbol 



Value 



total matter 

(in the units of critical density) 
non-relativistic matter 
baryonic matter 
relativistic matter 
vacuum-like energy 
Hubble constant 
radiation temperature 
helium mass fraction 
deuterium-to-hydrogen ratio 



A 
Y 



[D/H] 



^CDM + ^b 0.27 

0.045 

fi 7 + tt v 8.23 ■ 10" 5 
0.73 

70 km/s/Mpc 
2.725 K 
0.24 

1 3 • 10- 5 



1 



3 Hydrogen Feedbacks 

3.1 Semi-analytical Consideration of Resonant Feedbacks for Hy- 
drogen 

The n <-»• 1 transition rate is given by the following formula 



where A n [s _1 ] is the coefficient of spontaneous radiative n — > 1 transition, N n [cm -3 ] is 
the concentration of atoms in the state n, Ni [cm -3 ] is the concentration of atoms in the 
ground state, r]{v) is the number of photons per v mode, u n is the central frequency of 
n — > 1 line, g n , gi are the statistical weights of state n and the ground state. 

Use of Sobolev solution allows us to rewrite j2j) in more convenient form: 



where P n is the usual Sobolev escape probability, 77^" is the occupation number of radiation 
in the blue wing of n —> 1 line (see Fig. []]). 

Let us define the coefficient C n (the factor by which "ordinary" Sobolev transition 
rate (corresponding to the presence of proper resonant radiation only) is inhibited by the 
feedback effect): 



where J° is the n 1 transition rate in the case of equilibrium radiation in the blue 
wing of n <-> 1 line (i.e. 77+ = 77° where 77° is equilibrium (planckian) occupation number 
at frequency z/ n ). Neglecting induced n — > 1 transitions (in comparison with spontaneous 
n — > 1 transitions, since 77+ < 10~ 8 1 for any reasonable case) one can obtain: 






(3) 




(4) 



a 
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Nn ~ (gn/gl) ^1 

N n -(gJg 1 )r ] ^N 1 



(5) 
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Using quasistacionary solution of radiation transfer problem for n — y 1 line: 

V = f§(l-e-T)+V + e^ (6) 

and taking into account that for all Ly-lines optical depth r is much larger than unity at 
the central frequency and in the red wing, one can write: 

Vn^Vn - —vr (7) 

where rj~ is the occupation number of radiation in the red wing of n — > 1 line (see Fig. 
CD). Using (jTJ) one can obtain the following expression for C Q : 

C n = ^4 (8) 

Vn ~Vn 

In accordance with Zeldovich et al. (1968) and Peebles (1968) the populations of 
excited states are approximately determined relative to the state n = 2 by the Boltzmann 
distribution: 

N n g a ( E n -E 2 \ 

exp — — (9) 



N 2 92 V k vT 
Using (JTJ) and ([9]) it is easy to show (e.g. Peebles 1968) that: 



E n — E-) 



Vn = V2 6X P ( ^fr~ ) (10) 



or approximately 



la 

where r} a = rj 2 is the HI Lya occupation number (additional subscript a emphasizes the 
important role of HI Lya radiation for consideration of cosmological recombination). 

It is well known that in the absence of absorption the value r]+ can be found by 
using the following formula: 

^)=Ch(4) (12) 
where z' n for this case is given by the following relation 



z n = {l + z) — l = (l + z) — - _ n _ 2 1 (13) 



Using (TTT1) one can obtain from (PT2l) the following: 



r£(z) = rhM$Q (14) 



Substituting (fTTi ) and f fl~4l ) into (jHj) one can find that 



r _ (W^) - (^«)A£(Q) 
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Note that definition (fTBl results in rj® +1 (z' n ) = r]®(z). Using this one can obtain the 
following formula for C n : 

c " = mW^ (16) 

Let us define the following function: 

re = £-l (17) 

'la 

which has sense of relative overheating of Lya radiation in comparison with its equilibrium 
value (see top panel of Fig. [2|). Using this function the formula for C n can be rewritten 
in the following form: 

c n = i-r H «)/r H (z) ( 18 ) 



Formula ( 11811 shows us that in the first approximation the knowledge of the only 
function Tn(z) is enough for calculating all the feedback inhibition coefficients C n . 

Asymptotic of C Q at large n can be obtained by using Taylor's expansion for Tn(z' n ): 

r H «) ~ r H (z) + ^j^- « - z) (19) 

For the large n the following approximation is valid (e.g. Grin and Hirata 2009) 

z J n ~ z + 2(1 + z)n- 3 (20) 

Using (fT9l) and (l20l) one can obtain 

c 2(l + z)rflnr H _ 2 d\nT n 
" n 3 dz n 3 dln(l + z) 

Formula fl2~T]) is valid for principal quantum numbers n > 16 and redshifts z = 800 — 1800 
with relative accuracy within 10%. Dependence of dlnTn/d\n (1 + z) on z is shown in the 
bottom panel of Fig. [2l Note that this dependence is almost linear for z = 1000 — 1300. 
Taking this into account one can find simple approximation of Th for different values of the 
key cosmological parameters. Such an approximation may be useful for fast calculations 
of full set of inhibition coefficients C n without preliminary calculation of Th. 

Formula (!2~T1) shows that at any fixed z values C n decrease with increasing n by 
power law n~ 3 . This means that beginning from some value n max the total influence of 
transitions from high excited states with n > n max and feedback effect becomes negligible 
(at any given level of accuracy). 

It should be noted that C n are the functions of non-equilibrium radiation field at 
different moments z' (shifted relative current moment z). This radiation field depends on 
solution of cosmological recombination problem which in turn depends on C n in general 
case. This means that for maximal accuracy of numerical solution of considered problem 
one should use iteration method (Chluba and Sunyaev 2007, Switzer and Hirata 2008). 
Nevertheless it is evident that influence of taking into account C n on ionization history 
and non-equilibrium radiation field is small, i.e. correction to intensity of non-equilibrium 
radiation is much less than unity (no more than several percent). This allows us to use 
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unperturbed (i.e. calculated in the frame of standard TLA model without feedbacks) 
value of Lya radiation field for determination of function Th(z). 

Coefficients C n may be useful not only for improvement of TLA model but also for 
using in multilevel codes (without splitting on angular moment (Grachev and Dubrovich 
1991) and with splitting on angular moment (Burgin 2003, Rubino-Martin et al. 2006, 
Switzer and Hirata 2008, Grin and Hirata 2009)) for precise calculations of ionization 
history and recombination distortions of CMBR spectrum. As well as formalism of Sobolev 
escape probability, the formalism of inhibition coefficients C n allows us to avoid direct 
solution of radiation transfer equations for lines and transitions under consideration. 

Results of calculations of coefficients C n (z) are shown in the Figs [3] and EH 



3.2 Influence of Feedbacks on Hydrogen Recombination 

Knowledge of coefficients C n (z) allows us to calculate the feedback correction to the 
hydrogen ionization fraction Xmi and correspondingly to the concentration of free electrons 
AN e /N e . Ionization fraction xmi is determined by the standard TLA kinetic equation 
(Peebles 1968, Seager et al. 1999) with modified inhibition factor: 

n ^eflf + ^2sls / 00 n 

Phi + ^i e ff + ^2sis 

where A r eS is the new total effective coefficient of np<->ls transitions, A2 S i s is the coefficient 
of 2s— >ls two-photon spontaneous transition, /3hi is the effective total ionization coefficient 
from excited states of hydrogen atoms. The effective coefficient A r eS is given by the 
following formula: 



^ = E CnfPnA* exp -^-^ (23) 

where E n is the energy of level n. Taking into account the approximated expression for 
Sobolev escape probability 

g n A Q N m c 3 

one can find the correction to the total rate of np<-s-ls transitions in comparison with the 
"standard" 2p<-^ls rate: 



- »\ nT ^ (24) 



AJ 



2pls 



n>3 a \ a / 



(25) 



Including feedback 3^2 leads to the negative (i.e. decelerating) correction 
(1 — C2) while (n + 1) => n feedbacks (n > 3) lead to positive (i.e. accelerating) cor- 
rections. It has the following explanation: coefficients C n take into account not only 
re-absorption of (n + 1) — ► 1 photons in n — > 1 line (decelerating part of effect which is 
equal to {C n — 1)) but also n — > 1 transitions due to escape of photons from line pro- 
file (accelerating part which is equal to 1). Both these effects have not been taken into 
account in the standard TLA-model excepting 2p— »ls transitions due to escape of Lya 
photons from line profile. 
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Note that in traditional notation of recombination kinetic equation (e.g. Peebles 
1968, Wong et al. 2008) for allowance of hydrogen feedbacks one should replace usual 
Peebles -fTn-factor (i^H ld = c 3 / (87iu^H), where H is the Hubble constant as a function of 
z according formula (pp)) by the following: 

K*™ = K$ d (l + (26) 

V ^2pls / 

Results of calculations of relative change of free electron concentration during hy- 
drogen recombination epoch are shown in the Fig. El 

For comparison with result by Chluba and Sunyaev (2009) for hydrogen, the rel- 
ative difference of the free electron concentration AN e /N e between the Feedback TLA 
model and TLA model taking into account transitions from high excited np-states (n > 3) 
due to escape of Lyn photons from line profiles (e.g. Dubrovich and Grachev 2005) has 
been calculated (see Fig. [6]). The maximum of this relative difference is 0.259% at 
z ~ 1055, which is in very good accordance with corresponding result by Chluba and 
Sunyaev (2009). 



3.3 Effect of Neutral Deuterium 

The fact that deuterium Lya frequency is larger than hydrogen one means that redshifted 
hydrogen Ly/3 photons may be absorbed by neutral deuterium earlier than by hydrogen. 
This may lead to screening of redshifted hydrogen Ly/3 radiation from HI Is— >2p transi- 
tions by depth of neutral deuterium. Fung and Chluba (2009) have noted that existence 
of such screening depends on concrete conditions of radiation transfer, and estimate of 
screening effect may be very complicated problem in general case. In this subsection some 
necessary details of redshifted HI Ly/3 radiation transfer through DI Lya and HI Lya lines 
are considered. The absorption coefficients [cm -3 ] for hydrogen and deuterium in 2p — > Is 
line are given by the following formula 

Kx {u, z) = ^Ax,2p0X,2p(^)A r x,l s (27) 
9is 

where subscript X denotes component (X = H or D), 0x,2 P ( zy ) [Hz -1 ] is the absorption 
2p— >ls line profile for component X (J 0x,2 P {v) du — 1), Ax,i s is concentration of atoms 
X in ground state. 

Corresponding optical depths are given by the following formula 

T ^ z)= L s^ J wr d ^ (28) 

where z' = (1 + z)u' /u — 1. 

The value iV D)ls is controlled mainly by charge transfer reaction (Galli and Palla 
1998, Standi et al. 1998 and references therein): 

D + + H <-> D + H + (29) 

which tries to keep the fraction of neutral deuterium a little larger than the one of neutral 
hydrogen. This deviation between fractions of neutral hydrogen and neutral deuterium is 
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due to difference between their ionization energies which in turn is determined by isotopic 
shift for these atoms. This energy difference is ATdh = 43 K in temperature units. 
Therefore during hydrogen recombination epoch z = 800 — 1600 the relative deviation of 
neutral deuterium fraction from hydrogen one is an order about AT DH /T < 0.02 <C 1. 
Neglecting this deviation one can obtain the following formula (see also Wolf Savin 2002) 

iV D ,is^iV D ^ (30) 
iv H 

where iV D is the total concentration of deuterium atoms and ions, Nu is the same for the 
hydrogen. Using ([30]) one can estimate the absorption coefficient kd and optical depth 
td. 

Absorption coefficients and optical depths are shown in the Fig [7l The difference 
between central frequencies of absorption coefficients for hydrogen and deuterium is due 
to isotopic shift having relative value m e J (2m p ) ~ 27.2 • 10~ 5 . The difference between 
magnitudes of absorption coefficients is determined by the natural abundance of deuterium 
relative to hydrogen. The difference between thermal (Doppler) widths of absorption 
coefficients is due to different masses for hydrogen and deuterium atoms. 

From top panel of Fig. [7] one can see that at the central frequency of DI 2p— >ls 
line the absorption coefficient of deuterium is about an order of magnitude larger than 
hydrogen one (note that for hydrogen the wing of profile is at this frequency) . Neverthe- 
less, despite this one should take into account that regions and frequency ranges where 
the interaction between matter and radiation becomes significant are determined by the 
following condition r ~ 1. From bottom panel of Fig. [7] one can see that optical depth for 
deuterium achieves unity at the frequency where the optical depth for hydrogen already 
has the values much larger than unity (~ 10 4 ). Therefore, despite the fact that central 
absorption frequency for deuterium is larger than for hydrogen, the actual frequency of 
interaction for deuterium is less than corresponding one for hydrogen. This means that 
the redshifted Ly/3 photons of hydrogen interact with neutral hydrogen at transitions 
HI ls^2p earlier than with neutral deuterium at transitions DI ls^2p, i.e. effect of 
screening is absent in the considered case. 

4 Helium Feedback 

4.1 Semi-analytical Consideration of Feedback for Helium 

In the difference from radiation transfer during hydrogen recombination epoch when the 
medium is almost transparent for radiation at non-resonant frequencies (i.e. absorption 
is negligible) the transfer of redshifted Hel n— >1 photons during helium recombination 
epoch occurs at the presence of sufficient continuum absorption, namely absorption at 
ionization of neutral hydrogen atoms. Interaction of redshifted Hel n— >1 quanta with 
neutral hydrogen leads to the partial thermalization of this superequilibrium radiation. 
This should be taken into account in the solution of radiation transfer equation, so instead 
of formula (fl2l) one should use the following: 

Vn(z) = Cfl«) ex P (- r (n+l)^n) + Vn( z ) i 1 ~ e*P (-T(n+l)^n)) (31) 
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where T^ +1 \ =>n ^ s the optical depth due to ionization of neutral hydrogen. Since in present 
work we consider only 2 l P <-* l l S 2 3 P <-» l 1 ^ feedback for helium recombination we 
can write the following expression for "7"|ip^2 3 P : 

where z 2 3 P = (^ip/^p) (1 + z ) ~ 1 5 (V) is the ionization crossection of hydrogen 
ground state at frequency v' = i>2 3 p(1 + z')/{\ + z), Nm(z') is the concentration of 
neutral hydrogen atoms at the moment with redshift z' . Using Saha approximation for 
A/hi the optical depth t 2 I ip=>2 3 p can De estimated (with relative accuracy within 5% for 
z > 1800) by the following expression: 

H /x ca H (z^p) A/hi (z) ^b^ (2) / / In 7 h \\ /qq\ 

^p^p W ^ 777-^ t 1 " exp =- ^7— r - 7— tt-t (33) 



#(z) J H V V^B^(^23p) &b^(^ 

where In is the ionization energy of hydrogen ground state. 

Likewise the hydrogen case, let us define Thc = (^p/^p) — •"•) where ?7 2 ip is the 
actual occupation number in Hel 2 1 P — ► l 1 ^ line, i s the equilibrium occupation 
number in Hel 2 1 P —> l 1 ^ line. 

Substitution of (13T1) into (jSJ) gives us the following formula: 

n * THe(^2 3 p) / H \ /oa\ 

CW = 1 - r 7 s exp (-r 2 ip^ 2 3 P ) (34) 
Dependencies of C 2 3 P and tJi p ^ 2 3 P on redshift z are shown in the Fig. [3 



4.2 Influence of Feedback on Helium Recombination 

Knowledge of coefficient C 2 3 P (z) allows us to calculate the feedback correction to the 
helium ionization fraction XHeii = Nn e n/Nn e and correspondingly to the concentration of 
free electrons AN e /N e . Ionization fraction is determined by the standard TLA kinetic 
equation (Matsuda et al. 1969, Seager et al. 1999, Kholupenko et al. 2007, Wong et al. 
2008) with modified inhibition factor for ortho-channel of recombination: 

c _ (ggWggas) ^2 3 p (P 2 % + P&p) exp (-E 2 3 P2 3 S /k B T) 

Par + C 2 3 P (g&p/g&g) A 2*p (P^p + P 2 3p) exp (-P 2 3p 2 3s/ k B T) 

where A 2 zp is the coefficient of Hel 2 3 P — > l 1 ^ spontaneous transition, E 2 3 P2 3 S is the 
2 3 P — > 2 3 S transition energy, (3 or is the effective total ionization coefficient of helium 
ortho-states, P^p is the probability of Hel 2 3 P — > l 1 ^ transitions due to destruction of 
Hel 2 3 P — > l 1 ^ resonant photons by neutral hydrogen within line (Switzer and Hirata 
2008, Kholupenko et al. 2007, Rubino-Martin et al. 2008), P 2 r 3p is the modified Sobolev 
escape probability of Hel 2 3 P — ► l 1 ^ transitions (Kholupenko et al. 2008). Expressions 
for probabilities Pj! p and P 2 r 3p appropriate for TLA-models can be found in Wong et 
al. (2008) and Kholupenko et al. (2008). For comparison with Chluba and Sunyaev 
(2009), the dependence of the effective escape probability C 2 3 P (Pjp + P 2 3 P ) on redshift 
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z is presented in Fig. [9l This function is in very good accordance with result by Chluba 
and Sunyaev (2009). 

Dependencies of the free electron fraction normalized by total concentration of 
hydrogen N £ /Nh and correction AN e /N e on redshift z for the helium recombination 
epoch are presented in Figs [TO] and [EH correspondingly. 

In the Fig. [TT]for additional comparison with Chluba and Sunyaev (2009), AN e /N e 
in the case of partial taking into account hydrogen continuum absorption is presented 
("partial" means that radiation transfer within the resonant 2 1 P <-> l 1 ^ and 2 3 P <-> l 1 ^ 
lines has been calculated with taking into account the presence of neutral hydrogen, 
but radiation transfer between these lines has been calculated without the presence of 
neutral hydrogen). In this case the maximum of AN e /N e is 0.372% at z ~ 2018. Such 
perturbation of free electron number can affect the visibility function much larger than 
in the case of complete taking into account hydrogen continuum absorption (i.e. that 
would be important for analysis of damping tail of CMB anisotropy). Obtained result is 
in very good accordance with corresponding one by Chluba and Sunyaev (2009). Once 
again it should be noted that result for partial taking into account hydrogen continuum 
absorption is not valid and was obtained as intermediate one for comparison only. 

5 Results 

The main result of this paper is the semi-analytical approach (formulae [T7l IT8l l22l - l26 l [331 
- [35]) describing the feedback effect for resonant transitions in cosmological plasma. This 
approach allows us to calculate the free electron concentration N e taking into account 
feedback effect and correction AN e /N e relative to the "standard" free electron concentra- 
tion obtained in the frame of common TLA model (i.e. TLA model without taking into 
account transitions from high excited states n > 3 and feedbacks, e.g. recfast by Wong 
et al. 2008). 

The correction to the free electron concentration AN e /N e during hydrogen recom- 
bination epoch is presented in the Fig. [5l The maximal value of this correction is 0.216% 
at redshift about z ~ 1020. It can lead to the anisotropy power spectrum correction 
having the order about 0.1 - 0.2%. 

The correction to the free electron concentration AN e /N e during helium recombi- 
nation epoch is presented in the Fig. [TT1 The maximal value of this correction is 0.117% 
at redshift about z ~ 2308. Such value of correction is negligible at the current level of 
experimental accuracy (e.g. needed for treatment of Planck experimental data). 

All results obtained in this work (in the frame of modified TLA model) are in very 
good accordance with corresponding results obtained in the frame of multilevel model by 
Chluba and Sunyaev (2009). 

The main conclusions of present work are the following: 

1) The modification of ionization history of the Universe due to HI Ly(n + 1) =>- HI n«-»l 
feedbacks is important for correct analysis of CMB anisotropy power spectrum at the cur- 
rent level of experimental accuracy (e.g. Planck mission). This modification in particular 
can be taken into account by using semi-analytical approach developed in present paper. 

2) The modification of ionization history of the Universe due to Hel 2 1 P <-> l 1 ^ =>• 
2 3 P <-> l 1 ^ feedback is negligible for correct analysis of CMB anisotropy power spectrum 
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at the current level of experimental accuracy. The smallness of this modification is par- 
ticularly provided by screening of Hel 2 1 P <-> l 1 ^ resonant quanta from Hel 2 3 P <-> l 1 ^ 
transition by small amount of neutral hydrogen existing during helium recombination. 
It allows us to avoid accounting the feedbacks for helium at modeling of cosmological 
recombination. 
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Figure 1: Schematic picture of superequilibrium radiation transfer: occupation 
number r] as a function of frequency v is shown by thick solid line, absorption 
line profiles as functions of frequency v are shown by thick dashed line. Value 
v n is the central frequency of n — > 1 line profile, value u^ =1 is the frequency at 
which the optical depth in n — > 1 line achieves unity. 
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Figure 2: Top panel: Dependence of T H on redshift z. Bottom panel: 

Dependence of dlnT^/dln (1 + z) on redshift z. 
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Figure 3: Dependencies of C n on redshift z are shown by solid lines for differ- 
ent values of principal quantum number n indicated near lines. Dependence of 
(1 — C-i) on redshift z is shown by dashed line. 
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Figure 4: Dependencies of C n on principal quantum number n are shown by 
triangles (for z = 1500) and squares (for z = 1100). 
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Figure 5: Convergence of results at variation of number of levels involved in 
calculation. Curves presented are the relative differences between Feedback TLA 
model and common TLA model (without transitions from high excited states 
(n > 3) and feedbacks): result for n max = 2 (i.e. feedback 3 =>- 2 is taken into 
account) is shown by dashed line, result for n max = 3 (i.e. feedbacks 4 =>- 3 2 
are taken into account) is shown by dashed-dotted line, result for n max = 5 (i.e. 
feedbacks 6 =>- ... =>- 2 are taken into account) is shown by dashed-dotted line, 
result for n max = 30 (i.e. feedbacks 31 ... =>- 2 are taken into account) is 
shown by solid line. Note that curve for n max = 5 is very close to the curve for 

n max = 30. 
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Figure 6: The relative difference of free electron concentration between Feedback 
TLA model and TLA model taking into account transitions from high excited 
np-states n > 3 (i.e. TLA model with C n = 1). This result is obtained for the 
number of levels n max = 10. 
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Figure 7: Top panel: Is— >2p absorption coefficients for hydrogen (solid line) 
and deuterium (dashed line) as functions of relative frequency deviation from 
central frequency of hydrogen Lya-line. Curves are shown for moment z = 1300. 
Bottom panel: optical depths of absorption at HI Is— >2p transition (solid line) 
and at DI Is— >2p transition (dashed line). Dotted line shows the unity level. 
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Figure 8: Top panel: Coefficient C 2 3p as a function of redshift z with taking 
into account hydrogen continuum absorption (solid line) and without taking into 
account hydrogen continuum absorption (i.e. for 7" 2 ip=^2 3 P = ^' dashed line) 
Bottom panel: Optical depth T"^ P=>2 3 P as a function of redshift z (according 
formula (1331)). Dotted line shows the unity level. 
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Figure 9: Effective escape probability C 2 3p (P^ P + -P^p) as a function of redshift 
z with taking into account feedback (solid line) and without (i.e. for C 2 a P = 1, 
dashed line) 
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Figure 10: Dependence of N e /N H (relative number of free electron normalized by 
the total number of hydrogen atoms and ions) on redshift z: the result for model 
without feedback is shown by dashed line, the result for model with feedback and 
partial taking into account hydrogen continuum absorption (i.e. P|f p and P^ p 
take into account hydrogen continuum absorption, but t^ p=>2 3 P = 0) is shown 
by dashed-dotted line, the result for model with feedback and complete taking 
into account hydrogen continuum absorption (t^ip =>2 3 p ^ s calculated by formula 
(EH}) is shown by solid line. 21 




Figure 11: Relative change of free electron concentration AN e /N e as a function 
of redshift z: (i) the result for model with feedback and partial taking into 
account hydrogen continuum absorption (i.e. P^ p and Pjs p take into account 
hydrogen continuum absorption, but t^ p=>2 3 P = 0) is shown by dashed line, (ii) 
the result for model with feedback and complete taking into account hydrogen 
continuum absorption (t^i p=>2 3 p ^ s calculated by formula (133]) ) is shown by solid 
line. 
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